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Recently, the connection between Majorana fermions bound to the defects in arbitrary dimensions, and com¬ 
plex momentum roots of the vanishing determinant of the corresponding bulk Bogoliubov-de Gennes (BdG) 
Hamiltonian, has been established (EPL, 2015, 110, 67005). Based on this understanding, a formula has been 
proposed to count the number (n) of the zero energy Majorana bound states, which is related to the topological 
phase of the system. In this paper, we provide a proof of the counting formula and we apply this formula to a 
variety of 1 d and 2d models belonging to the classes BDI, Dill and D. We show that we can successfully chart out 
the topological phase diagrams. Studying these examples also enables us to explicitly observe the correspon¬ 
dence between these complex momentum solutions in the Fourier space, and the localized Majorana fermion 
wavefunctions in the position space. Finally, we corroborate the fact that for systems with a chiral symme¬ 
try, these solutions are the so-called "exceptional points", where two or more eigenvalues of the complexified 
Hamiltonian coalesce. 
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1. Introduction 

Topological superconductors Q] are systems which can provide the condensed matter version of Ma¬ 
jorana fermions, because they can host topologically protected zero energy states at a defect or edge, 
for which the creation operator (j) :=0 ) is equivalent to the annihilation operator (je= o). These localized 
zero-energy states obey non-Abelian braiding statistics (2] [3], which can find potential applications in 
designing fault-tolerant topological quantum computers (2]0). Although Majorana fermion bound states 
have not yet been conclusively found in nature, they have been theoretically shown to exist in low di¬ 
mensional spinless p- wave superconducting systems (2j |5j, as well as other systems involving various 
heterostructures with proximity-induced superconductivity which are topologically similar to them (6f 

ED- 

Non-interacting Hamiltonians for gapped topological insulators and topological superconductors, in 
arbitrary spatial dimensions, can be classified into ten topological symmetry classes Il3j|l5l . character¬ 
ized by certain topological invariants. Moreover, there exists a unified framework for classifying topo¬ 
logical defects in insulators and superconductors (16), which follows from the bulk-boundary correspon¬ 
dence and identification of the protected gapless fermion excitations with topological invariants char¬ 
acterizing the defect. Here we focus on Id and 2d Bogoliubov-de Gennes (BdG) Hamiltonians with the 
particle-hole symmetry (PHS) operator squaring to +1, which can be categorized (TD into three classes: 
BDI, Dill andD. 

In our earlier work (l7), we have explored the connection between the complex momentum solu¬ 
tions of the determinant of a bulk BdG Hamiltonian ( A/bcK,) in arbitrary dimensions, and the Majorana 
fermion wavefunctions in the position space associated with a defect or edge. We have found that the 
imaginary parts of these momenta are related to the exponential decay of the wavefunctions, localized 
at the defects, and hence their sign-change at a topological phase transition point signals the appearance 
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or disappearance of Majorana zero mode(s). Based on this understanding, we have proposed a formula 
to count the number ( n ) of the zero energy Majorana bound states, which is related to the topological 
phase of the system. This formula serves as an alternative to the familiar Z and Z 2 topological invariants 
GUSH EGO and other counting schemes I19H22I . 

In this paper, we prove this formula and apply it to a variety of Id and 2d models belonging to the 
classes BDI, Dill and D. We show that we can successfully chart out the topological phase diagrams. Study¬ 
ing these examples also enables us to explicitly observe the correspondence between these complex mo¬ 
mentum solutions in the Fourier space, and the localized Majorana fermion wavefunctions in the position 
space. Finally, we also corroborate the fact that for systems with a chiral symmetry, these solutions can 
be identified with the so-called “exceptional points” (EP’s) I23H29I . where two or more eigenvalues of the 
complexified Hamiltonian coalesce. EP’s are singular points at which the norm of at least one eigenvector 
vanishes, when certain real parameters appearing in the Hamiltonian are continued to complex values, 
and the complexified Hamiltonian becomes non-diagonalizable. The concept of EP’s is similar to that of 
a degeneracy point, but with the important difference that all the energy eigenvectors cannot be made 
orthogonal to each other. In previous works, EP’s have been used l36tl'3~4l l to describe topological phases 
of matter for Id topological superconductors/superfluids. 

The paper is organized as follows: in section[2] we review the results obtained earlier (17) for counting 
the number (n) of Majorana zero modes bound to defects, based on the bulk-edge correspondence. In 
section[3] we provide a proof of the counting formula. In section[4] we consider some Id and 2d models 
in the class BDI and apply the EP formalism to count n. Section[5]is devoted to the study of edge states for 
Hamiltonians in class Dill, where we illustrate the applicability of EP solutions as the chiral symmetry 
exists. In section [6] we discuss some systems in the class D and conclude that EP’s cannot be related 
to the Majorana fermion wavefunctions for such Hamiltonians, because chiral symmetry is broken. We 
conclude with a summary and outlook in section[7] In appendix[A] we provide a simple example to show 
how one should choose the correct EP solutions such that their imaginary parts are continuous functions 
in the parameter space in order to evaluate our counting formula. 


2. Counting formula for the Majorana zero modes 


In this section, we review the connection GD between the complex momentum solutions of 
det[H B dG(k)] = 0, and the Majorana fermion wavefunctions in the position space associated with a de¬ 
fect or edge. 

We consider a topological defect embedded in (or at the boundary of) a d-dimensional topologi¬ 
cal superconductor. Let m be the dimensions of the defect, parametrized by the Cartesian coordinates 
r_L = iri,...,r d -m) and ry = (r d _ m+ i,...,r ri ), located at rj_ = 0. Let k ± = k ± Cl = (fci,..., k d - m ) and k y = 
[k d - m+ \,...,k d ) be the corresponding conjugate momenta, where k± = |kxl and Cl is the unit vector 
when written in spherical coordinates. 

For a generic Ff B dG, l et k A and k i O' - !»■■■> Q) be the two sets of complex k ±-solutions for 
detjiTsdGtk)] = 0, related by jlm(fc^)} = -{Im(fcg)}, after k± has been analytically continued to the complex 
plane. One should be careful to choose solutions such that their imaginary parts are continuous functions 
of the parameter(s) which tune(s) through the transition, and the solutions in one set are related to the 
other by changing the sign of their imaginary parts throughout. This point has been illustrated by an 
example in appendix [ a| Assuming the Majorana wavefunction to be of the form ~ exp(-z[rj_|) in the 
bulk, the correspondence iA;j_ «— -z has been established (l7). At a topological phase transition point, one 
or more of the Im(fc^ /B )’s go through zero. When Im(fc^ /B ) changes sign at a topological phase transition 
point, the position space wavefunction of the corresponding Majorana fermion changes from exponen¬ 
tially decaying to exponentially diverging or vice versa. If the former happens, the Majorana fermion 
ceases to exist. A new Majorana zero mode appears in the latter case. The count ( n ) for the Majorana 
fermions for a defect is captured by the function 


/({Aihky.A) 


1 

2 


E ( si S n { Im [ k m (Wzh k n , Cl) ] } - sign jlm [ fc{ /B ({A°}, k{j, n°)]} ] 

1=11 ' 


( 2 . 1 ) 
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where ({A,},k||,O) are the parameters appearing in the expressions for fc^ /B , and ({A9},k^,0°) are their 
values at any point in the non-topological phase. 

If there is a chiral symmetry operator & which anticommutes with the Hamiltonian, the latter takes 
the form 

^chiralC k ) = ( o 5 )> (2.2) 

in the momentum space, for the corresponding bulk system with no defect. On analytically continuing 
the magnitude k^ = k- |k| to the complex kj_ -plane, at least one of the eigenvectors of flchirafik) collapses 
to zero norm where 

det[.e/(k)] = 0 or det [,s/^(k)] = 0. (2.3) 

These points are associated with the solutions of EP’s for complex k ±-values where two or more energy 
levels coalesce. Furthermore, these coalescing eigenvalues have zero magnitude since det[^(k)] = 0 (or 
detts/^k)] = 0) also implies det[Tf C hi ra i(k)] = 0. Hchiraitk) becomes non-diagonalizable, as in the com¬ 
plex k± -plane, det[.s/(k)] = 0 =e> det^^fk)] = 0 (or vice versa). However, at the physical phase transition 
points, the imaginary parts of one or more solutions vanish, and det[srf (k)] = detU^fk)] = 0 for those so¬ 
lutions, making H C hir a i(k) once again diagonalizable and marking the disappearance of the corresponding 
EP’s. 

Since it satisfies equation j2.3) , each EP solution corresponds to a Majorana fermion of a definite 
chirality with respect to 0. If ^(k) = x^ T (-k) holds, then the two sets of EP’s are related by {k ] A } = -{fc B }, 
one set corresponding to the solutions obtained from one of the two off-diagonal blocks. In such cases, 
the pairs of the Majorana fermion wavefunctions are of opposite chiralities. 


3. Derivation of the counting formula 

A simple derivation of the counting formula in equation 113 can be motivated as follows: 

1. Let us consider one of the solutions given by j = 1. In the non-topological phase, say phase “0”, 
fci(U°},k{|,n c ') gives no Majorana zero mode and hence does not give rise to any decaying mode 
localized at a defect. On the other hand, in a topological phase, say phase “t”, with a Majorana 
wavefunction ~ exp[-|Im(A^)| r±], fc^({A(},k||,Q)| phaset localized at r± = 0 and zero at r± = oo, 
should now give rise to an admissible decaying zero mode solution. This implies that there is a 
change in sign of Im(fc^) from -1 to +1 when we jump from phase “0” to phase “t”. 

2. Majorana zero modes must occur in pairs, though they might be localized far apart. Hence, if 
biphase c corresponds to a Majorana mode localized at rj_ = 0, then A: B | phase t must correspond to 
one localized at r_L = oo, where fc B = (k A ) *. Hence, whether or not we are in the topological phase “t” 
is captured by the function f\ = \ [sign{lm[fc^ B ({A ! -},k||,n)]}-sign{lm[fc^ B ({A9},kpfl 0 )]}| taking 
the value 1 or zero. 

3. From the above discussion, it may seem that the counting formula should be given by 

\ E | :si S n { Im [■ k m k ll--■*)]} - si S n { Im [■ k m ( {A ° k °’ ]} | ■ 

However, this is not quite correct. To understand this, let us consider the scenario when at least two 
of the solutions, say k A and k A are such that sign[Im(A:^)] | phase0 = -sign[Im(fc^)] | phase0 - This implies 

that in the trivial phase, the wavefunction given by c\ ex P(i A^| phase0 rj + C 2 ex P(i^ilph a seo rj -) 
inadmissible for not being capable of satisfying the boundary conditions — the only solution is 
ci = c 2 = 0. In another topological phase, say “t”, let sign[Im(fc^)] | phase{ = -sign[Im(fc^)] | phase0 and 
sign[Im(fc2)]| phasef = -sign[Im(fc2)]| phase0 . This means that both Im(fc^) and Im(fc^) change sign 
when we jump from phase “0” to phase “t”. However, they still should not give any Majorana 
zero mode in the phase “t”, because Cj exp(i A^| phase f r±) + c 2 ex P(i^Alph a set ?J -) cannot satisfy the 
boundary conditions. So, the correct formula is given by equation s 


33703-3 




I. Mandal 


4. EP formalism for the BDI class 

In this section, we consider some Id and 2d spinless models in the BDI class, which can support 
multiple Majorana fermions at any end of an open chain. For systems in this class, there exists a chiral 
symmetry operator ©, such that H B dG can be rotated to the form Ff C hirai in equation (2.2) . 

After reviewing the transfer matrix scheme to find Majorana fermion solutions localized at an edge, 
we show how EP solutions in the complex kj_ -plane can be used to count the number of Majorana zero 
modes in a given topological phase. We also make emphasis on the connection of these EP solutions with 
the position space wavefunctions calculated in the real space lattice with open ends. 


4.1. Transfer matrix approach 

Kitaev (2]| suggested the model of a Id p- wave superconducting chain, which can support Majorana 
zero modes at the two ends. For a finite and open chain with N sites, the Hamiltonian takes the form 




(4.1) 


where p is the chemical potential, w and A are the nearest-neighbour hopping amplitude and supercon¬ 
ducting gap, respectively. The pair of fermionic annihilation and creation operators, Cj and cj, describe 

the lattice site j, and obey the usual anticommutation relations {c/, cj} = 0 and {Cj, cj,} = Sjji . The Majo¬ 
rana mode structure of the wire can be better understood by rewriting the above Hamiltonian in terms 
of the Majorana operators 


a j — c j + Cj, bj = -i jet - Cj j, 

satisfying 

aj-a\, bj-b\, {aj,bj'}- 0, {aj,aji} - {bj,bjj — 28 jji . 

Then, the Hamiltonian reduces to 


(4.2) 


j tv j tv-i 

H K=-~Lv a jb ] --Y 1 [iw-A)ajbj+i-{w + A)bjaj+i]. (4.3) 

Z 7=1 Z 7=1 

This chain can support one Majorana bound state (MBS) at an edge for appropriate values of the 
parameters. More recently, a variation of the model was considered with next-nearest-neighbour hop¬ 
ping and pairing amplitudes (35). A general version of such longer-ranged interactions with all possible 
hoppings and pairings was studied I36II37I with the Hamiltonian 

j tV q N-q 

Hi — — — pcij bj — i ^ ^ [J-r a j bj+ r + Jr dj+r bj ], (4.4) 

4 7=1 r= 1 7=1 


where the /+ r ’s are real parameters, and 0 < q < N. These models can support multiple MBSs at an edge. 
If we impose periodic boundary conditions (PBC’s), the Hamiltonian can be diagonalized by a Bogoliubov 
transformation: 


H i = ~L( c l c ~ k ) h i {k) ( c\)’ 

I V ( Jrcos(kr) -i/ r sin(fcr) \ 

1 r k q { iJrsm(kr) -J r cos (kr) j’ 



(4.5) 


where the anticommuting fermion operators (cj, Cfc) are suitable linear combinations in the momentum 

space of the original (c/, ct) fermion operators. The energy eigenvalues are given by 

J J 


Edk) = ± 2 


12 r ,2 

2j/ r cos(fcr) + 2j/ r sin(A:r) . 


(4.6) 
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We now review the transfer matrix approach l35j)38l to identify the number of MBSs at each end of the 
chain for this model. The transfer matrix can be obtained from the Heisenberg equations of motion for 
the Majorana operators in equation l |4.4) : 

. d«; . 4 . d bj 4 

2i——= -i £ J- r bj +r , 2\——-'i £ J r a,j +r . (4.7) 

m r=-q r=-g 

Assuming the time-dependence to be of the form aj - Aj e~ lE ‘ r and bj - Bj e~ lElt , the £/ = 0 (zero energy 
modes) are given by the recursion relation of the amplitudes: 

q q 

£ J-rb j+r = 0, £ J r a ]+r = 0. (4.8) 

r=-q r=-q 

Clearly, it will suffice to solve one set of the recursive equations to obtain the solutions for both. Assuming 
Aj - A ] a and Bj - Ag, we get the polynomial equations 

£/ r Af r = 0, £ /_ r Af r = 0. (4.9) 

r=-q r=-g 

An MBS can exist if we have a normalizable solution, i.e., if |A A | < 1 or |A B | < 1, if the solution is to be 
localized at the left end. Similarly, for a mode to be localized at the right-hand end of the chain, we must 
have |AaI > 1 or |A B [ > 1. Depending on the number of constraint equations (or boundary conditions on 
the amplitudes), one should determine the number of independent MBSs at each end of the chain. 


4.2. Relation of the EP formalism with the transfer matrix approach 

Let us apply the EP formalism Ifl7l[34l to the Hamiltonian in equation |4.5| . First we rotate it to the 
off-diagonal form 


hi l0d m = u]hi{k)Ui = y B “ [k] 


Alik) 




q q 

Aiik) ——2 £ [ Jr cos {kr) +ij r sin(fcr) ], B[{k)--2 £ [j r cos(fcr) - i/ r sin(fcr) ]. (4.10) 

r=—q r=-q 


The EP’s where either Ai[k) or _B/(fc) vanishes, are given by the solutions 

4 

£ } r A£ f =0, where 

r—-q 
4 

£ J- r Ag£ = 0, where 

r=-q 

Comparing equations ( |4.9| l, (4.11| and |4.12| , it is easy to see that the solutions for EP’s in the complex 
fc-plane for the PBC’s correspond to the MBS solutions for the open boundary conditions (OBC’s). Since 

I^aZ/b/I < 1 Im(A: A ;/ B z) > 0 o |Aa/bI^1> (4.13) 

IAaz/bzI > 1 => Im(fc A // B z) < 0 o |A a /b I > 1 > (4.14) 


A a z = exp (ifc A /), 

(4.11) 

Abz = exp (ifc Bi ). 

(4.12) 


a sign change of Im(/c A // B /) indicates a topological phase transition, by which we move from a phase 
where an MBS can exist to the one where that particular zero mode gets destroyed. This is related to 
the fact that Im(fc A // B z) , s are related to the exponential decay of the MBS position space wavefunctions 
localized at one end of the open chain. 

Choosing Jo = -j, ]\ - J2- J-i = J-2 = and all other J r ’s to be zero, we can get a system 
supporting up to four Majorana zero modes at each end of the chain. The phase diagram obtained using 
equation is shown in figure [T] 
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Figure 1. (Color online) The topological phase diagram of the Hamiltonian described by equation (4.4) , 
with Jq = -j, Ji = J 2 — yy, J-i = 7-2 = ^ 2 ^, and all other / r ’s set to zero. Here, n labels the number 
of Majorana zero modes at each end of the chain, as captured by the function /(p, A) defined in equa¬ 
tion (2.1) . 


Instead, for the parameters Jo — -y cos(02)> h = cos(0i), /_ 1 = sin((/>i), and all other TV's set 

to zero, we get a system having three EP’s for either Ai[k) = 0 or B[{k) - 0. For this model, up to two Ma¬ 
jorana zero modes can appear at an edge. The phase diagrams for JIM — 0.625 and JIM - 1.3, obtained 
using equation o are shown in figure [2] 

We should note another important point: if there are Q EP solutions for either Ai{k) - 0 or Bi{k) — 0, 
clearly there are 2 Q solutions in total. However, for counting the zero modes in equation (2.1) , we should 
consider only one set, where the two sets obey the relation 

^-Al — 1/A.bz or ^ai — —koi. (4.15) 



J/M=1.3 

liiiii 1 1 1 1 1 



0 1 2 3 4 5 6 


01 

(b) 


Figure 2. (Color onl ine) Panels (a) and (b) show the topological phase diagram of the Hamiltonian de¬ 
scribed by equation |4.4|, with Jo = -ycost^l, /1 = cos(0i), /-1 = sin(</>i), and all other J r ’s set 

to zero. Here, n labels the number of Majorana zero modes at each end of the chain, as captured by the 
function /(p, A) defined in equation H3- 
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As we have already seen, these two sets correspond to the wavefunctions of the MBSs at the two opposite 
ends. Evidently, the MBSs exist in pairs at the two ends and the topological phase is characterized by their 
number at each individual end. 


4.3. Single-channel ferromagnetic nanowire 


The Id Hamiltonian for a ferromagnetic nanowire embedded on Pb superconductor (39) with a single 
spatial channel (i.e., no transverse hopping) is given by 

Hm = XXhjvti*)'*'*, ¥* = (c k] ,c kl ,c f _ k ) T , 

k 

— £(k)&oTz + [A s ob + ApSin(fc)d- a] t x + V- ar Tq , £(A) — —2t cos(fc) - p. (4.16) 


Here, k is the Id crystal momentum, 'Pfc is the four-component Nambu spinor defined in the particle- 
hole (t) and spin (tr) spaces, and V is the Zeeman field which can be induced by ferromagnetism. Also, 
A s and A p are proximity-induced s-wave and p -wave superconducting pairing potentials, respectively, 
with d determining the relative magnitudes of the components of the p -wave superconducting order 
parameter A a p ( a,p = t,|). In our calculations, we use d = (1,0,0) and V = (0,0,1/). This Hamiltonian 
belongs to the BDI class with the chiral symmetry operator given by 0 — a x T y . 

The eigenvalues of the Hamiltonian are given by: 


Ei (fc) = ± yj (fc) + V 2 + Af + A p sin 2 (k) - e\ , E 2 {k) — ± \J$ 2 (k) + V 2 + A z s + ApSin 2 (fc) + e\ , 

e\ - 2 J V 2 [A 2 + f 2 (fc)] + A 2 Ap sin 2 (A;) . 


(4.17) 


A level crossing can occur if either Ei(k) = 0 or E 2 (k) - 0. However, for a finite V and A s , the latter is 
impossible. Hence, a level crossing takes place when E\ (fc) = 0 for k = 0 or n for the appropriate values 
of the parameters, which also indicates that this corresponds to the appearance of zero energy modes. A 
generic complex value of k corresponding to E\ ( k ) = 0 can be obtained by solving 

[y 2 -£ 2 (fc) +Ap S in 2 (A;)] 2 + 4£ 2 (A;) A 2 sin 2 (fc) = 0, V= \]v 2 -A 2 s . (4.18) 

We can rotate the Hamiltonian in equation (4.16) to the chiral basis, where it takes the form 


H^\{k) = Ul h N] (fc) U 2 




f -1 - i 

0 

1 + i O' 

( o \ 

i h?\k) 0 )’ 

1 

%= 2 

0 

0 

1 + i 

1-i 

0 -1-i 

0 1-i 



l 1-i 

0 

1-i 0 , 


h'T Ik) = 


-£(fc) + iA p sin(fc) - V A s 

A s (A:) - iAp sin(A) - V )’ 


h? ] (k) 


-£{k) -iApSin(fc) - V A s 

A s £(fc) + iAp sin(fc) - V )' 


(4.19) 


If either detth^ 1 (fc)] = 0 or det [h^ 1 (A;)] = 0 for a complex fc-value, this leads to the vanishing of the norm 
of one of the four eigenvectors of H c ^(k), signalling the existence of an EP for that value of k. At an EP, 
H^ik) is thus non-diagonalizable. 

The solutions for the EP’s are given by either 


or 


det[/z„ (A:)] = 0 => t/ 2 -^ 2 (fc) + A 2 sin 2 (A;) = — 2i<f (A:) A p sin(A;) 

=> k= k“ S2 = -ilnj siV- p + s 2 \/(.si V- p) 2 -At 2 (2f+Ap) _1 


det [ Az, 1 (A:)] = 0 


V 2 - £ 2 {k) + A 2 sin 2 (A) = 2i^(fc) Apsin(fc) 


k = *4,s 2 = “ilnj SiV-p + s 2 \J [si V-pf-At 2 {2t-A p ) 


(4.20) 


(4.21) 
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Figure 3. (Color online) Parameters: A s = A p = 0.1 f, V = 1.5 1 corresponding to the Hamiltonian in 
equation |4.16} . (a) Energy bands Eiglk), given in equation {4.17| , have been plotted in blue and red, 
respectively, as functions of p/ 1 . (b) Plots of Im (l“ lS2 ) versus p/ 1. (c) /(p) giving the count of the chiral 
Majorana zero modes as a function of /i/ 1. 


where (si = ±1, S 2 = ±1). Clearly, k = fc“^ 2 also solves equation i4.18|, which corresponds to two coincid¬ 
ing zero energy solutions (where two levels coalesce (27) for a complex /c-value). 

The plots of the energy bands, Im(fc“ S2 ), and /(p) have been shown in figure [ 3 ], using the values 
A S -Ap- O.lf and V - 1.5 1. 

Now, let us try to understand the existence of the EP’s throughout a given topological phase and 
their disappearance right at the phase tran sition points, the latter being tied to the sign change of the 
Im(fc" / j )’s. The Hamiltonian in equation 4.19 , when written in position space, gives the following 


equations for the Majorana zero modes, yr + = (u+,0) T and 
spectively): 


= (0, U-Y (with chirality +1 and -1, re- 


d\ + /i + A p d x - V 
A, 

d 2 x + / 1 - A p d x -V 


A s -d 2 x -[i-\ p d x -V 


A s 


-d x - fi + A p d x - V 


u+ = 0, 


U— — 0 . 


(4.22) 


Here, we have assumed a continuum for an open wire and set t — 1. For i//_, let us assume the trial 


solution U- = ^exp {-z r x) 

r 

deJ Z r + V-* P Zr-V 


u' 


U r 


. The complex z r ’s must satisfy the quartic equation 


-zr r - fi+ Ap z r - V 


= 0 


[Zr+H-ApZr] =V , 


(4.23) 
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whereas for small k, from equation j4.21| , we get 

[-k 2 + p + iA p k) 2 = V 2 , (4.24) 

indicating correspondence ik -z r . The magnitude of z r will determine the admissible MBS solutions 
subject to OBC’s (as analyzed in an earlier work @0}), just as in the transfer matrix analysis for the Id spin¬ 
less lattice case. Hence, here also we have been able to establish the relation between the existence of EP’s 
in the complex fc-plane (for the periodic Hamiltonian) and the localized Majorana zero modes at the ends 
of an open chain. 


4.4. Two-channel time-reversal-symmetric nanowire system 

MBSs in a two-channel TRS nanowire proximity-coupled to an 5-wave superconductor have been re¬ 
cently studied OED- The low-energy model for the lowest bands of the system is described by the effective 
Id 4 x 4 BdG Hamiltonian: 

H N 2 = 2>I/»iV2 M'Pfc. n = (C fct ,C tl ,ct -cl ) T , 

k 

k 2 

h N2 = Zik)a 0 T z + v(ka z -p c a 0 )T x + B-(TT 0 , £(fc) = -- p, (4.25) 

2m 

where p c is the momentum when the gap closes, B is a magnetic field for the Zeeman term, and ( v, p ) are 
effective parameters. We have set p c = 2 vm for our calculations. Since this nanowire system belongs to 
the BDI class when B is perpendicular to the spin-orbit-coupling direction, we will take B = {B, 0,0) in our 
analysis. Then, the chiral symmetry operator is given by & — n z T y . 

The eigenvalues of the Hamiltonian are given by: 


E\ (7c) - ±\j B 2 + k 2 v 2 +4 m 2 v 4 + £ 2 (fc) - e 2 , E 2 {k) = ± y B 2 + k 2 v 2 + 4 m 2 v 4 + ^ 2 (fc) + e 2 , 


e 2 = 2 y4 m 2 v 4 ( B 2 + k 2 v 2 ) + B 2 £ 2 (fc). (4.26) 

We can have two levels coalescing if E\ (fc) = 0 for a complex A:-value obtained by solving 

[B 2 + k 2 v 2 - ^ 2 (fc) - 4 m 2 v 2 } 2 + 41 2 (A:) k 2 v 2 = 0. (4.27) 

As before, we rotate the Hamiltonian in equation |4.25f to the chiral basis, where it takes the form 


*0*1 = h N2 {k ) U 3 = 


o K 2 {k) 

h^ 2 {k) 0 


1 

lh = - 
2 


I 0 -1 - i 0 1 + i 

0 1 - i 0 1 - i 

1+i 0 -1-i 0 

1-i 0 1-i 0 


h N2 {k) = ( -hk)+i[kv + 2mv 2 ) 


-l{k)-\[kv + 2mv 2 ) 


B 


-f(fc) +i[kv-2mv 2 ) 


hj[k) = 


B 


B 


-0,{k) -\[kv-2mv 2 ) )' 


The solutions for the EP’s are then given by either 

det \hu 2 {k)\ = 0 => [l{k) - ikv^ - B 2 - Am 2 v A 


k= fc“ 1>S2 = —iln si \j m 2 v 2 - 2mp+ 2\s 2 m VAm 2 v 4 - B 2 + i mv\, 


or 


det [h^ 2 {k)\ — 0 


[l{k) + ikv] 2 = B 2 -Am 2 v 


„2„4 


k - k\ l S2 = -iln|si \jm 2 v 2 -2mfi + 2is 2 m \/~Am 2 v 4 - B^- imv^, 


(4.28) 


(4.29) 


(4.30) 
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£i, 2 (k) ^1,2 (k) 




- 3 - 2-10 1 2 3 


P/(v p c ) 

(c) 


Figure 4. (Color online) Parameters: v = 1, m = H{2v 2 ), p c = corresponding to the Hamiltonian in 
equation |4.25| . Panels (a) and (b) show the energy bands £ 12 (fc), given in equation j4.26| , as functions 
of p, for B = 0 and B = 3, respectively. £ 1,2 (fc) have been plotted in blue and red, respectively. Panel (c) 
shows the contourplot of f{p) giving the count “n” of the MBSs in the pl[v p c ) - Bl{vp c ) plane. 


where (si = +1, s 2 = ±1)- Clearly, k = also solves equation 1 4.27' and hence corresponds to the coa¬ 
lescing of two energy levels at the zero value in the complex /c-plane. Choosing v = 1 and m = 1/(2 v 2 ), 
the energy bands for B = 0 and B = 3, and the contourplot for f[p,B) [defined in equation b have 
been shown in figure [ 4 ] Once again we find that f{p,B) gives the correct topological phase diagram in 
figure[4](c). Needless to add that here also exp(i fc^^l’s determine the admissible solutions for the MBS 
wavefunctions in the position space, at the ends of an open chain. 


4.5. Majorana edge modes for the Kitaev honeycomb model 

In this subsection, we consider the EP-formalism for a 2d lattice Hamiltonian in the class BDI. The 
Kitaev honeycomb model f42l can be mapped onto free spinless fermions with p -wave pairing on a hon¬ 
eycomb lattice, using the Jordan-Wigner transformation. The solutions for the edge modes for a semi¬ 
infinite lattic^] have been studied earlier I49H52I . The momentum space Hamiltonian in terms of the 

1 We would like to point out that this system is different from two-dimensional p x +ip y fermionic superfluids, whose excitation 
spectra include gapless Majorana-Weyl fermions (43). Volovik showed that in such chiral superfluids, the fermionic zero modes 
along the domain wall have the same origin as the fermion zero modes appearing in the spectrum of the Caroli-de Gennes-Matricon 
bound states in a vortex core SID- This correspondence can be understood by picturing the chiral fermions as orbiting around 
the vortex axis, analogous to the motion along a closed domain boundary. 
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Majorana operators is 


H H = L(% ^W k )(g), h/,M = ( k = ik Xl ky) 


A(k) = -2i 

(k x k y ' 
h + Ji cos 1 ^ 

+ J 2 cos 

,kx + ky jj 

+ 2 

. ( k x k y ' 
,lSm ( 2 . 

+ J 2 sin 

,k ’ +k, )\ 

B (k) = 2i 

1 k x — k y 
J 3 + h cos 2 

+ J 2 cos 

kx + ky j 

+ 2 

. ( k X - ky 

2 J 

+ J 2 sin 

kx+ ky j 


with the eigenvalues 


(4.31) 


E{ k) = +2 


h + h cos 


kr k\ 


+ h cos 


k x + ky 


Ji sin 


k x ky 


+ h sin 


kx + ky 


1/2 


.(4.32) 


We will consider two kinds of edges l49]l50l . namely, zigzag and armchair, which can support Majo¬ 
rana fermions. We will find the phase diagram using the EP’s corresponding to these edges setting either 
A(k) = 0 or 5(k) = 0, after complexifying the momentum component perpendicular to the edge. For this 
2d case, / in equation is a function of U\,J2,h, fc||), where k\\ is the momentum along the Id edge 
being considered. 

One can have a zigzag edge in the y-direction, according to the convention of Nakada et al. (49], so 
that we will complexify k± = k x , and k\\ — k y will be one of the parameters determining the topological 
phase transition points. The solution for B{k x - k±, k y = k\\) = 0 is given by 


fc_i_ = -2iln 


Ji exp(ifc||/2) +J 2 exp(-ifc||/2) 

h 


(4.33) 


f(J]— Jl -(3, ky) 




f ( Ji +J2 ky) 

f(Ml — ^21>*^3 1 ky) 


k x ) 


fUi=Ji=Ji or J t <J 2 , k x ) 

lr 


0 

(c) 


2 


2 



0 

(d) 


2 


2 


ji 


k x 


Figure 5. (Color online) Topological phase diagrams for edges for the 2d honeycomb lattice described by 
equation |4.31 , as captured by the function /(/i,/ 2 ,/ 3 , k\\) defined in equation E3- Panels (a) and (b) 
show the number of chiral Majorana zero modes for a zigzag edge, while panels (c) and (d) show the same 
for an armchair edge located at the top of a semi-infinite lattice. 
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Majorana zero modes exist for all the values of lc y if J\ + J 2 < h -There is no edge state if |/i - J 2 \ > h- For 
Ji= J 2 = J 3 , edge states exist if \k y \ > 2tt/ 3. These results have been plotted in figures[5](a) and[5](b). 

For the armchair edge (49) in the x-direction on the top of the lattice, we will complexify k± - k y , and 
fc|l = k x will be now one of the parameters determining the topological phase transition points. The two 
EP’s for A{k x - fc||, k y = k±) - 0 are given by 


k y = -2iln 


~h exp(-ifc||/2)± yj / 2 exp(-ifc||) -4 /i/ 2 
2/2 


(4.34) 


No Majorana zero mode exists for any value of k x if J] - I 2 = h or Ji < h- For J\ > I 2 , a Majorana 
fermion can exist for a specific range of values for k x . Figures [ 5 ] (c) and[s](d) show these topological 
phases, obtained using equation (13- 

Equations |4.33) and (4.34| are seen to coincide with the solutions of the Majorana edge states ob¬ 
tained earlier by the transfer matrix formalism 14911521 . 


5. EP formalism for the Dill class 

A point defect in class Dill can support a Majorana Kramers pair (MKP) corresponding to doubly de¬ 
generate Majorana zero modes, whereas a line defect can support a pair of helical Majorana edge states. 
Both are characterized by a Z 2 topological invariant. The chiral symmetry operator © can be defined 
such that the Hamiltonian in class Dill can be brought to the block off-diagonal form [equation (13], just 
like for the class BDI. 


5.1. Id model 

A simple Id model of topological superconductivity in the class Dill is described by the Hamiltonian 

ED 


H m i = Y J *[. hmi m'Vk, n = (c fcT , c k j ,cl kl - c f _ k] ) T , 

k 

h m i ik) = [^ m i(fc)CT 0 + ARsin(fc)o- z ]T z + A cos(k) cj 0 t x , £ m i(fc) = tcos(fc)-p. (5.1) 

This system may be realized in a Rashba wire that is proximity-coupled to a nodeless s+ wave supercon¬ 
ductor. The energy eigenvalues are given by: 


E\(k) = + - A r sin(fc)] 2 + A 2 cos 2 (A;), E 2 [k) = ± \J[£, m \{k) + A R sin(fc)] 2 + A 2 cos 2 (A;), (5.2) 

whose plots are shown in figure[6](a) for A = 0.1 1 and A R = 2f, as fil t is varied along the horizontal axis. 

Observing that a chiral symmetry operator © = cro^y exists in the presence of Ji z , we rotate the 
Hamiltonian in equation (53 to the chiral basis, where it takes the form 




f 0 

-i 

0 

i 

<j(fc) = u\ hmi {k) U 4 = ( h £ [k) h ’» 1' [ ' C) ) , 

£ 

II 

Sth 

0 

-i 

1 

0 

0 

i 

1 

0 



l 1 

0 

1 

0 


h“ n 1 (Ac) = diag|iAcos(fc) -^,„i(fc) + A R sin(fc), iAcos(fc) -£ m i(fc) - A R sin(fc)j, 
h\ nl (k) = diag|- iAcos(fc) - ^ m i(fc) + A R sin(fc), -iAcos(fc) - f m i(fc) - A R sin(fc)|. (5.3) 


The solutions for the EP’s are then given by either 

det[h ; “ j(fc)] = 0 => iAcos(A:) - Z m i(k) - siA R sin(A;) 

H+s 2 y^u 2 - A|- (t — i A) 2 


k=k Z% = ~ iln 


t - i (A + si A r ) 


(5.4) 
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£i, 2 (k) 



(a) 

Figure 6. (Color online) Parameter: A = O.lf corresponding to the Hamiltonian in equation H3- (a) En¬ 
ergy bands Ei^lk), given in equation {5.2} , have been plotted in blue and red, respectively as functions 
of p/ f, for Ar = 2 1. (b) /(p/ f, Ar/ t ) giving the count “n” of the chiral Majorana fermions. 



or 


det | h ml {k)\ =0 


k=k^l=-\\n 


iAcos(fc) + Zmiik) = SiA R sin(fc) 


M+ *2 \Jn 2 -A,^-{t + iA) 2 

t + i (A + Si Ar) 


(5.5) 


where (si = +1, S 2 = ±1). For this model, we note that the t wo diffe rent sets of EP solutions, related by 
{Im(fc)}| set=A = -{Im(fc)}| set=B , are obtained from equations 1 5.4' and' 5.5 when we set si = 1 and si = -1, 
respectively. This is related to the fact that [h“ x (A:)]"*" ^ [fc“ x (— A:)] T (where h l ml [k ) = [h“ n 1 (fc)] t ) for real k. 
However, we have argued before that for equation (23 to work, we must take all the EP solutions from 


one of the sets related by a negative sign of Im(fc). Using either (fc^f“ , ) or (fc^y“ , k™f l ) (rather than 

both), figure [ 6 ](b) gives the correct topological phase diagram in the p/t- A R /f plane, for A = O.lf. We 
clearly see that there exist phases with a pair of MBSs, which correspond to one Kramers doublet (MKP). 


5.2. 2d model 

The Id model of a Rashba semiconductor combined with a nodeless s+ wave superconductor can be 
easily generalized to a 2d system, described by the Hamiltonian (54) 

H m2 = IX h m2 W = (Ckf. Cfcj, , -C^ ) T , 

k 

hm 2 (k) = [<fm 2 (k) + A m (k)] CT 0 T Z + 2A R [sin(fcJ a y - sin(fc y ) <j x ]t z , 

£ m2 00 = -2t[cos(fcJ + cos(fcy)] - p, A m (k) = A 0 + 2 A 1 [cos(k x ] + cos(fc y )], (5.6) 

where A m (k) is the s± wave singlet pairing potential that switches its sign between the centre (0,0) and 
the corner (n,n) of the 2d Brillouin zone, when 0 < |Ao| < 4Ai. 

The energy eigenvalues are given by: 


E\ (k) 

E 2 ( k) 


= ± ]j $ m 2 (k) - 2 Ar \J sin 2 (A : x ) + sin 2 (fc y ) 
= + \j £ m2 (k) + 2 A r y/sin 2 (k x ) + sin 2 (fc y ) 


+ A?„(k), 


+ Am(k), 


(5.7) 


whose plots are shown in figure[7](a) for A R = A 0 = Ax = 2t, as yu /1 is varied along the horizontal axis. 


33703-13 
























I. Mandal 



Figure 7. (Color online) Panel (a) shows the plot of energy levels E) i2 (k) of equation H3 in blue and red, 
respectively, for the Hamiltonian in equation (5.6) , as functions of /i/ 1. We have used Ar = Ao = Ai = 2f. 
Panels (b) and (c) show the contourplots of /(p/f, Ao/Ai, k y , Ar = 2 1, Aj = t) in the Ao/Ai - p/f plane, 
giving the count “n” of the Majorana edge states along the y-direction, for fc y = 0 and k y = n, respectively. 


Rotating the Hamiltonian in equation (5.6) to the diagonal basis of the chiral symmetry operator 6 = 
a 0 T y , we get 


H%W = ulh mM U, = ( h ,° K 

u iA m (k)-f m2 (k) 2A R [sin(fc y )-isin(fc x )] \ 

mA 1 1 2AR[sin(fc y )+isin(fc*)] iA m (k) - f m2 (k) )’ 

h l ,,, = ( -iA m (k)-£ m2 (k) 2A R [sin(fc y )-isin(fc x )] \ 

m2 { 2A R [sin(fc y ) + isin(fc^)] -iA m (k) -< m2 (k) )' 

The equations for the EP’s, corresponding to edge modes along the y-direction (so that k\\ = k y and k± = 
k x ), for k y — 0 and k y = n, are given by 

det [h“ 2 ( k )] I, =° => {iA m (k)-^ m2 (k)}| fc (Q , =2siA R sin(fcJ, (5.9) 

I ICy —(Uy/TJ / 

and 

det [ h l m2 (k) 1 I = 0 => {iA m (k)+£ m2 (k)}| fc (07t) = 2siA R sin(fc x ). (5.10) 

L J I Ky— \\JfTt) y 
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The solutions for k y - 0 are: 


2f + jU + i(2Ai+A 0 ) «2 yW - A 0 ) (4 Ai + A 0 - - 4i t) - 4 A^ 

2 1 + 2 i (Ai + Ar) 2 f + 2i (Aj + Ar) 


and 


-2r-/i + i(2Ai + A 0 ) \/-(i// + A 0 )(4Ai + Ao + ip + 4ir)-4A| 

---+ — - - 

2 1 — 2 i (A^ + Ar) 2 1 — 2 i (A^ + Ar) 


Those for k y -n are: 


and 


2r-^i + i(2A 1 -A 0 ) $2 yf (A 0 -i/t) (4A x - A 0 + i/J-4ir) -4A^ 

2 £ + 2i (Ai + Ar) 2f + 2i(Ai + S\ Ar) 


2t— u + i(A 0 - 2 Ai) «2 J (Ao + i^) (4 A x - A 0 - i/j + 4i t) - 4 A^ 

- 11 - + —- - 

2 f — 2 i (A^ + Ar) 2 1 — 2i (A^ + Si Ar) 


(5.11) 


(5.12) 


(5.13) 


(5.14) 


Here, fc"' 2 5 “ and correspond to the vanishing of det[h^ 2 (k)] and det[/jJ );9 (k)], respectively, and 
(si = ±1,S2 = ±1). Two distinct sets of EP solutions, related by {Im(fcjJ}| set=A = -{Im(fcx)}| set=B , are ob¬ 
tained by setting s\ = 1 and s\ = -1, respectively. Using either (£” 2 “, k+\ l S2 ) or (fc™, 2 “, fc™ 1 ^) (rather 
than both) in equation j2.1) , figures [7] (b) and[7](c) give the desired topological phase diagrams in the 
Hit- Ao/Ai plane, for Ar = 21 and Ai = t. The topological phases with a pair of helical Majorana edge 
states are clearly seen. 


6. Broken time reversal symmetry: class D 

In this section, we consider Id and 2d Hamiltonians in the symmetry class D, where the TRS is broken. 
We will see that the EP formalism in the complex /cj_ -plane is not applicable for such systems. 


6.1. Id spinless model 

We examine the spinless model described by the Hamiltonian 


N -1 , 

h d\ = E 

i=i 


- w 


c] Cj +1 + A Cj Cj+i - w* ct +1 Cj + A* Cj +1 Cy] - £ H (cj Cj - , 


(6.1) 


which looks similar to Hk in equation |4Hj, but with the important difference that the TRS is broken by 
the fact that w and A can be complex numbers (37). Without any loss of generality, we can choose A to 
be real and encode the entire phase-difference (0) between A and wq by writing w — wq e 1 ^, where wq is 
real and positive. 

With PBC’s, one can write the corresponding BdG Hamiltonian in the momentum space as: 


#di = -E( c l c -k ) h m (k)l j ), 

k V L - k I 

hoilk) = [2 M/osin(0) sin(fc)] To - [2 wqcos(([> ) cos(fc) + /i] t z + 2 Asin(fc) T y . (6.2) 

The energy eigenvalues are given by: 

E(k) = 2 Wq sin((/>) sin(fc) ± e t , e t - y[ 2 Wq cos i<p) cos(fc) + /i]“ +4A 2 sin 2 (fc). (6.3) 
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Hence, it follows that two levels become degenerate when 

St- 0 => 2 wo cos(t/>) cos(fc) + p = + 2iAsin(fc). (6.4) 

When extended to the complex k- space, there can be EP’s where the norm of an eigenvector van¬ 
ishes. For convenience, we rotate hoiik) to find the points (EP’s) where the Hamiltonian becomes non- 
diagonalizable: 


h m , od(fc) = C/ f 


D 1 


hmik) Udi = 


2 wo sin(0) sin(fc) Ad i (fc) 

Boiik) 2 wo sin(t/>) sin(fc) 

1 -1 


Udi - —p i , , 

v/2 l - 1 - 1 

Ajji(k) -2 Wq cos(0) cos(fc) + p + 2iAsin(fc), 

Boiik) -2 wq cos(</>) cos(fc) + p - 2iAsin(fc). 


(6.5) 


The EP’s are given by (A:) = 0 or (k) = 0. At such points, two levels coalesce for a complex value of 
k satisfying equation j6.4) . 

However, we immediately observe that these EP’s do not correspond to zero energy modes, which 
appear when 

det[/iDi(fc)] = 0 => [2 wq cos(0) cos(fc) + yu] 2 + 4A 2 sin 2 (fc) = 4 Wq sin 2 ((p) sin 2 (fc) (6.6) 


<p =0 



/'/W'o 

(a) 


<p=2n/5 



- 3 - 2-10 1 2 3 


///*v 0 

(c) 


0=7T/1O 



///w 0 

Cb) 


0=JT/ 2 



- 3 - 2-10 1 2 3 


///W 0 

(d) 


Figure 8. (Color online) Panels (a), (b), (c) and (d) show the contourplots of /(p, A) corresponding to 
the Hamiltonian in equation (6.2} , giving the count “n” of the MBSs in the p/tt'o - A/ wq plane for <p = 
0,7r/10, 27r/5 and^/2, respectively. The blue regions have n = 1 MBS at each end of the open chain, while 
the purple regions correspond to the n = 0 trivial phases. 
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is satisfied. We note that since det[h m (fc)] is equal to the product of the energy eigenvalues, vanishing of 
det[hoi(fc)] implies the condition for the existence of a zero energy solution. 

Although the EP description no longer applies now to the existence of MBSs, we can find the complex 
fc-values satisfying 

det [hDi,od(fc)] = 0 => 2 wo cos (</>) cos (k) + p = ±2isin(fc) \J . A 2 - Wg sin 2 ((p). (6.7) 

We can solve for the fc-values either with the “+” or the sign on the RHS (rather than both), and plug 
in the roots of that equation into the formula in equation The function /(p, A) will still give the 
number of MBSs in a given topological phase. Once again we emphasize that to count the zero modes 
in equation s we should include only one of the two sets of roots related by a sign change of Im(fc), 
as these two sets correspond to the wavefunctions of the pair of MBSs at the two opposite edges. We 
have shown the plots of /(p, A) in figure[8]for four different values of (f>. It indeed captures the correct 
Z 2 topological invariant (n = 0, or 1). No zero mode exists, i.e., the system is entirely gapped in certain 
regions in the (f> - p plane where Im (fc) vanishes, as Im (/c) is related to the exponential part of the MBS 
wavefunction in the real space. 


6.2. 2d spinless model 

The following Hamiltonian gives a model of a p + ip wave superconductor on a square lattice f55l : 


Hd2 = ~Z( 4 C -k Ul*(k)( ^ ]. 

k v ' V C_ k / 

h D2 (k) = [2 t x cos (k x ) + 2ty cos(fcy) - p] t z + d x sin(fcj t x + d y sin(fc y ) r y , (6.8) 

where [t x , t y ) are the hopping strengths and (d x ,d y ) are the pairing amplitudes along the (x,y)-direc- 
tions, and p is the chemical potential. The energy eigenvalues, given by 


E{ k) 


= ± \J [2 t x cos(k x ) + 2 t y cos(fcy) - p] 2 + dlsm 2 {k x ) + d y sin 2 (fc y ), 


(6.9) 


are plotted in figure[9](a) for t x - t y - d x - d y - 1, as p is varied along the horizontal axis. 

We consider the edges parallel to the y-axis, so that k\\ = k y and k± = k x . Complexifying k x , we 
can compute the number of non-chiral Majorana fermions propagating along these edges with momenta 
k y - 0 and k y = n, using equation d2.1|. The complex k x -values satisfying det[h B2 (k)] = 0, for k y = 0 and 
k y = jt, are given by 


2 t x cos(k x ) + 2 t y - p = i s 1 d x sin {k x ) 

• 2 5l (1 - ty) + s 2 ^4(fy- l) 2 + d x -4 t x 


k x = k SliS2 =-iln 


d x + 2 Si t x 


( 6 . 10 ) 


and 


2 t x cos(k x ) - 2 t y - p = i si d x sinfe) 

" 2 (1 + ty) + S 2 \J 4 [ty + 1)^ + d x — 4 t x 


k x = k SliS2 = —i In 


d x + 2 Si t x 


( 6 . 11 ) 


respectively. Here, (si = ±1,S2 = ±1), and we need to use either si = 1 or si = -1 (rather than both) in 
equation to obtain the phase diagrams shown in figures[9](b) and[9](c). 


7. Conclusion 

We have established the relation of the EP solutions for complexified momenta to the Majorana 
fermion wavefunctions bound to a topological defect in a system with a chiral symmetry, by studying 
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Figure 9. (Color online) Parameters: d x = d y = 1 corresponding to the Hamiltonian in equation 6.8 . 
Panel (a) shows the plot of energy levels E{k), given in equation |6.9| , as functions of p, for t x = t y = 1. 
Panels (b) and (c) show the contourplots of f[t x , ty, k y , p = 2) in the t x - t y plane, giving the count “n” of 
the non-chiral Majorana edge states along the y-direction, for k Y = 0 and k y = n, respectively. 


some explicit examples in Id and 2d. These models include both spinless and spinful cases. We have 
shown that such EP solutions cannot exist for systems in class D, where there is no chiral symmetry. The 
generic formula, which was proposed earlier (T7) to count the number of Majorana zero modes in arbi¬ 
trary dimensions, has been demonstrated to chart out the desired topological phase diagrams for the wide 
variety of systems we have considered. The detailed study of these models also helps us illustrate how 
one distinct set of complex /cj_-solutions for det[H B dG(k)] = 0, related by {Im(/cj_)}| set=A = -{Im(fcJ}| set=B , 
should be used while using our formula. An explicit proof of the counting formula has also been dis¬ 
cussed. For a system with or without a chiral symmetry, the imaginary parts of these solutions in the 
complexified k j_-plane are related to the exponential decay of the Majorana fermion wavefunctions in 
the bulk in the position space. Hence, the imaginary parts of h of the solutions in one set undergoes a 
change of sign across a topological phase transition point, if the number of Majorana zero modes at a 
defect changes by n. 
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-4-2 2 4 
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Figure 10. (Color online) For edges of the system described by equation A.l , corresponding 
(a) Plots of versus p. (b) Plots of Im(ifcj s ,) versus p. (c) /(p) giving the count of the 

zero modes as a function of p. 


to ky = 0: 

Majorana 


A. Choice of EP solutions 

In this appendix, we provide a simple example to show how one should choose the correct EP solu¬ 
tions such that their imaginary parts are continuous functions in the parameter space for our counting 
formul;(5 Let us take the 2d class D Hamiltonian (56): 

H m = Y,[ 4 c- k ) ^D3(k) ^ , /to 3 (k) = k x T x + k y T y + (k 2 -h)t z , (A.l) 

k v ' V c _k J 

where /i is the chemical potential. The system is known to be topological for p > 0 and non-topological for 
p < 0. For non-chiral Majorana fermions along the edges parallel to the y-axis with momentum k y - 0, 
one should solve for det[ho 3 Uc x , k y = 0)] = -k x - [k x - p) 2 = 0. We will have four solutions which can be 
written as either 


or 


k 


nc 

5l,5 2 


= S 1 



\A~4 M 
2 


+ hi 


= (kj-n)+isik x = 0 


kx - k c Sl ,s 2 


-isi + s 2 i/4/i- 1 
2 


(A.2) 


(A.3) 


where 


det[h m {k x ,k y = 0)] = 


(A.4) 


and (si = +1, S 2 = ±1). 

2 We thank Victor Gurarie for suggesting to clarify this point. 
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I. Mandal 


The plots of ImlA:”^,), ImCfc^ S2 ) and f(fi) have been shown in figure[lO| We find that Im[fc” c _(ju)] = 
-Im[fc” e _(/i)] and Im[fc" c + (/r)] = -Im [£”_+(//)]. Also, Im[fc^ + (ju)] = -Im[A;£ _(//)] and Im[fc^ _(/u)] = 
-Im[fc£ + Qu)]. So naively, one might think that from the solution set in equation dA.2b, we can evalu¬ 
ate /(/i) by using any one of the four pairs given by (A:”JL, k™ + ), {k__, k_^ + ), {k”°_,k™ + ). 

Whereas, from the solution set in equation i |A.3) , f{[i) is expected to be obtained by using any one of the 
four pairs given by (fc£ + , k+_), {k c + + , fc£ + ), (fc£ _, k c + _), ( k c _ _, k c _ + ). One can check that this is true for the 
set in equation ^A.3) , rather than for the one in equation jA.2| . In other words, one gets the wrong phase 
diagram on using the first set. This is because all the Im(fc”^ S2 )’s are not continuous functions of p. 
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niflpaxyHOK 3B'fl3am/ix CTaniB Maftopam/i 3 Bi/iKopi/icTaHHflM 
KOMn/iexcHi/ix iMny/ibciB 

I. MaHfla/i 

iHCTHTy TeopeTHHHoi cf)i3HKi/i 'TlepuMeTp", Baiep/ioo, OHTapio N2L 2Y5, KaHafla 

HeipoflaBHO (EPL, 2015, 110, 67005) 6yao BCTaHOB/ieHO 3B'n30x Mi>x (jtepMiOHaMU ManopaHH, 3B'n3aHHMi/i 3 
fle^exTaMH y flOBiabHixi BUMipHOCTi, i KOMn/ieKCHUMM iMnynbCHnMn xopeHnMM fleTepMiHaHTa BiflnoBiflHoro 
o6'EMHoro raMinbTOHiaHy 6oro/uo6oBa-fle XeHa. Ea3yK)HHcb Ha pbOMy po3yMiHHi, 3anponoHOBaHO t()opMy/iy 
flan niflpaxyHKy Hi/icaa (n) 3B'n3aHwx CTaHiB MawopaHH 3 HyabOBOK) eHepriero, nxi noB'n3am 3 TonoaoriHHOK) 
t|)a30K) cttCTeMH. B pin CTani flaeTbcn BMBifl tf>opMyan niflpaxyHKy, nxa 3acTOCOByeTbcn flo hh3kh Id i 2d mo- 
fleaeti, uto HaaexaTb flo KaaciB BDI, Dill i D. riOKa3aHO, nx MO>KHa ycniLUHO no6yflyBaTH TonoaoriHHi cf>a30Bi 
fliarpaMH. BHBneHHn flaHHX npHKaafliB fl03Boane nBHO cnocTepiraTH BiflnoBiflHicTb Mi>K phmh KOMnaexcHHMH 
p03B'n3KaMH flan iMnyabcy b cpyp'e npocTopi i a0Kaai30BaHHMH XBHabOBHMH (JjyHKLtinMH <|)epMioHiB MatiopaHH 
b no3HLtiiiHOMy npocTopi. HaKiHepb, niflTBepflxeHO tfiaKT, tpo flan CHCTeM 3 xipaabHOto CHMeTpieio pi po3B'n3KH 
e Tax 3B3HHMH "BHHnTKOBHMH TOHxaMH", fle flBa hh 6iabuje BaacHHX 3HaneHb ycxaaflHeHoro raMiabTOHiaHa 
3aHBaK)TbCn. 

KmoHOBi cnoBa: eMHarxos/ tohki/i, (pepMiOHM Mafiopam/i, BDI, Dill, D, niflpaxyHOK 
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